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THEORETICAL STUDY OF THE TRANSONIC LIFT OF A DOUBLE-WEDGE PROFILE WITH 

DETACHED BOW WAVE ^ 

By Walter G. Vincenti and Cleo B. Wagoner 


SUMMARY 

A theoretical study is described of the aerodynamic character- 
istics at small angle of attack of a thin, double-wedge profile in 
the range of supersonic flight speed in which the bow ware is 
detached. The analysis is carried out within the framework of 
the transonic (nonlinear) small-disturbance theory, amd the 
effects of angle of attack are regarded as a smaR perturbation on 
the flow previously calculated at zero angle. The mixed flow 
about the front half of the profile is calculated by relaxation 
solution of a suitably defined boundary-value problem for the 
transonic small-disturbance equation in the hodograph plane 
(i. e., the Tricomi equation) . The purely supersonic flow about 
the rear half is found by an extension of the usual numerical 
method of characteristics. Analytical results are also obtained, 
within the framework of the same theory, for the range of speed 
in which the bow wave is attached and the flow is completely 
supersonic. 

The calculations provide, for vanishingly small angle of 
attack, the following information as a function of the transonic 
similarity parameter: (1) chordwise lift distribution, (S) lift- 
curve slope, and (S) position of center of lift. As in previous 
studies, the aerodynamic characteristics of a profile of given 
thickness ratio show little variation with free-stream Mach num- 
ber as the Mach number passes through 1. As the Mach 
number is increased to higher values, however, the lift-curve 
slope rises to a pronounced maximum in the vicinity of shock 
attachment and then declines. Correspondingly, the center of 
lift moves forward toward the leading edge and then returns aft. 
These findings are in marked contrast to the behavior of the drag 
coeflicient at zero angle of attack, which was found in earlier 
work to decrease monotonicaUy as the Mach number increased 
above 1 . At Mach numbers above that for shock attachment, the 
results of the present calculations tend toward those given by 
classical linear theory. 


INTRODUCTION 

The theoretical problem of the transonic flow over a thin, 
double-wedge profile at zero angle of attack has been treated 
in several papers iu recent years. These papers have in com- 
mon that they employ the simplifying concepts of the 
transonic smaU-disturhance theoiy and utilize the hodograph 
transformation to liueaxize the resulting mathematical prob- 
lem. FoUowmg this approach, Guderley and Toshihara 
(ref. 1) began by solving the problem for a free-stream Mach 
number of 1, using analytical methods for the mixed flow 
over the front wedge and the method of characteristics for 
the purely supersonic flow over the rear. Somewhat later, 
the present authors, using a combination of relaxation 
methods and the method of characteristics (ref. 2), extended 
the results to free-stream Mach numbers greater than 1, 
where a detached bow wave occins ahead of the profile. At 
about the same time. Cole (ref. 3) obtained an analytical 
solution for the flow over the front wedge at subsonic flight 
speeds, utilizing, in effect, the special assumption of a vertical 
sonic line from the shoulder of the wedge. More recently, 
Trilling (ref. 4) has been able to remove this special assump- 
tion and, with the aid of less stringent approximations re- 
garding the flow over the rear wedge, to extend the solution 
for the subsonic case to include the complete profile. As a 
result of these investigations, the problem of the double- 
wedge profile at zero angle of attack may be regarded as 
substantially solved within the limitations of the transonic 
snqaU-disturbance theory. The experimental studies of 
Liepmann and Bryson (refs. 5 and 6) and Griffith (ref. 7) 
indicate that the theoretical findings are in fundamental 
agreement with the physical facts. 

In a recent paper (ref. 8), Guderley and Toshihara have 
continued their investigations of the double-wedge profile at 
Mach nmnber 1 by considering the influence of a vanishingly 
small angle of attack. The basic idea in this later work is to 
regard the effects of angle of attack as a first-order perturba- 
tion on the nonlinear flow previously calculated at zero angle. 
This approach leads to a linear boundary-value problem in 
both the physical and hodograph planes. The calculation for 


1 8nperso<leaNAOATN283a,“TlieoretIcaIStudyofttieTian*onloLUtofaI>ouMs-W0dgBP)rofllewltliDetao!iedBowWave”by Walter a. Vlncentland Oleo B. Wagoner, 1052. Portions 
ol this work were also reported at tlia vmth International Congress on Theoretical and Applied Mechanics, Istanbul, Turkey, August 2(1-28, 1952, 
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the front wedge is still carried out, however, in the hodograph 
plane, since the basic procedures can then be taken over 
directly from the previous work. By this means, Guderley 
and Yoshihara obtain results for the lift-curve slope of the 
profile at zero angle of attack and for the corresponding 
distribution of lift along the chord. 

The aim of the present paper is to extend the results for 
the double wedge at angle of attack to the case of supersonic 
flight mth detached bow wave. The fundamental ideas of 
Guderley and Yoshihara are followed in reducmg the calcu- 
lations for the front wedge to a perturbation problem in the 
hodograpb plane. The detailed formulation of the problem 
is, however, necessarily different in the present case. The 
boundary conditions for the problem appear in terms of the 
results already obtained at zero lift (ref. 2), and the solution 
is carried out by numerical methods which differ only slightly 
from those devised for the earlier work. The lift on the rear 
wedge is calculated by an extension of the method of charac- 
teristics. The body of the paper is devoted to the detailed 
formulation of the boundary-value problem in the hodograph 
plane and to a discussion of the fiffal results. Noteworthy 
differences between the numerical procedures used in the 
present work and those already described in reference 2 are 
treated in appendices at the end of the report. 
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NOTATION 


PRIMABY SYMBOLS 


critical speed (i. e., speed at which the speed of 

flow and the speed of sound are equal) 

numerical coefficient 

(See eqs. (39) and (40).) 

airfoil chord ^ 

Tj.. /e • i lift per unit span 

lift coefficient, — 

2-c 


moment coefficient for moments taken about 

, ,. , moment per unit span 

Ic&ding ^ 


pressure coefficient, ^ 

integral defined by equation (45) 
numerical constant (See eq. (10).) 

Mach number 

slope of segment of Mach line in characteristics 
net 

static pressure 

local lifting pressure (i. e., difference between 
static pressures on upper and lower surfaces) 
dynamic pressure 
airfoil thickness 
speed of flow 
Cartesian coordinates 
generalized Cartesian coordinates 
(See eq. (43).) 


chordwise position of center of lift 


slope of curve of lift coefficient versus true 
angle of attack evaluated at zero angle 

slope of curve of moment coefficient versus true 
angle of attack evaluated at zero angle 


a normalized angle of attack; also denotes true 

angle of attack when used in derivatives 

such as etc. 
da 

/3 absolute value of at left-hand limit of lattice 

y ratio of specific heats (7/6 for air) 

A basic lattice interval 

C fimction of ij and 6 

(See eq. (A6).) 

ij, 5 normalized speed of flow 

(See eq. (la).) 
rii, rji special values of 

(See fig. 20.) 

6,6 normalized inclination of flow; 6 also denotes 

true inclination of flow in equation (lb) 

(See eq. (lb).) 

6u normalized half-angle of wedge 

transonic similarity parameter 
(See eq. (13).) 
p fluid density 

\f/ stream fimction 

6^'^ iucremental values of stream function 

(See eqs. (A9) and (All).) 

SUBSCRIPTS 

points in characteristics net 
(See fig. 19.) 

components of total stream function 
(See eq. (39).) 
copditions in free stream 
singular solution 

(See eqs. (A6) and (A7).) 
value at a prescribed lattice point 
conditions at critical speed 

SUPERSCRIPTS 

( ) quantities determined at zero angle of attack 

( )' derivative with respect to normalized angle of 

attack evaluated at zero angle 

BOUNDARY-VALUE PROBLEM IN HODOGRAPH PLANE 
PESCHIPnON OF FLOW FIELD 

Pigure 1 is a drawing of the idealized, inviscid flow which 
may be expected about a wedge profile when the angle of 
attack is sufficiently less than the semiapex angle of the wedge. 
Pigure 2 shows the corresponding hodograph representation 
'of the flow over the front wedge, which is the region of prime 
theoretical concern. Except for the substitution of the 
detached bow wave in place of the infinite free stream, these 
representations follow the lines assumed by Guderley and 
Yoshihara in reference 8. The corresponding draivings for 
zero angle of attack, which are fundamental to the present 
case, have been described in detail in reference 2. 

In the present example, the path of the central streamline 
in the physical and hodograph planes is briefly as follows: 
The streamline leaves the bow wave in the physical plane 
(or the shock polar in the hodograph plane) at point A. It 
then proceeds with decreasing subsonic speed to a stagnation 
point O on the imderside of the profile. At 0 the streamline 
branches. The lower branch runs downstream along the 
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Figure 1. — Flow about 'double-wedge profile at angle of attack. 

lower surface of the profile with fixed iaclination but in- 
creasing speed. The sonic speed is reached at the shoulder 
L, where the speed then increases discontinuously in accord 
with the Prandtl-Meyer relations. The shoulder itself maps 
in the hodograph onto the upgoing characteristic LM. 
The upper branch of the central streamline proceeds from O 
upstream along the surface of the profile. The inclination 
here is again fixed by that of the surface, and the speed 
increases to the sonic value at the leading edge J. At this 
point the flow is characterized by another Prandtl-Meyer 
expansion to supersonic speed. 


Figure 2. — Flow about front wedge in hodograpb plane. 

The flow configuration which should be assumed on the 
upper surface near the leading edge is open to conjecture. 
Since the geometrically available angle of turn will, for any 
thin airfoil, be greater than the 130° perpiissible for expansion 
to a vacuum, a region of separation is to be expected. If the 
angle of attack is not too great, this region will probably 
bo closed, with the central streamline reattaching to the 
upper surface a small distance behind the leading edge. 
This reattachment wUl be followed by a compression of the 
flow through a system of shock waves whose arrangement is 
sketched only formally in the physical plane (and not at all 
in the hodograph plane, where the correct representation 
would probably he on several sheets). The effects of the 
flow near the leading edge will be mentioned later, but the 
exact process wUl remain undefined. Whatever the details, 
the speed on the upper surface will return to a subsonic 
value at some point K just downstream of a terminating, 


normal shook wave. IVom K the central streamline con- 
tinues at fixed inclination downstream along the upper 
surface, the speed increasing once more to the sonic value 
at the shoulder B. At this point another expansion takes 
place, similar to that which occurs at the corresponding 
point on the lower surface. In this case the shoulder is 
represented in the hodograph by the downgomg character- 
istic BG. 

The supersonic expansion fan from the shoulder at B (and 
similarly at L) is discussed in detail in reference 2. Suffice 
it here to say that the supersonic flow field, of which the 
expansion fan is the initial part, is separated into two regions 
by the Mach line GE, which runs from the shoulder to the 
sonic point on the how wave. (This line was termed the 
“separating” Mach line in reference 2.) The supersonic 
flow in the region upstream of the Mach fine GE is interde- 
pendent with the subsonic field between the how wave and 
sonic line. To obtain a solution for the front wedge, a prob- 
lem in transonic flow must therefore be solved for the 
subsonic field and the interdependent portion of the super- 
sonic expansion fan. Conditions in the supersonic flow 
downstream of the Mach line GE have no influence upon 
the subsonic field. The continuation of the flow beyond GE 
can be accomplished by purely supersonic methods once the 
solution of the transonic problem is known. 

Aside from the obvious lack of symmetry in the present 
case, the main difference between the flow here and that 
previously studied at zero angle of attack is the existence in 
the present problem of the localized supersonic region in the 
vicinity of the nose. As pointed out, conditions in this 
region are difficult to formulate. The problem has been 
considered by Guderley and Toshihara (ref. 8) in the course 
of their work at Mach number 1. They find that, if the nose 
region is disregarded in the hodograph and the boundary 
condition along BIB is fulfilled all the way in to 0, then the 
influence on the lift of the resulting flctitious flow at the 
nose is of somewhat higher than the second order in the 
angle of attack. This suggests that the effects of the real 
flow at the nose may be neglected in a flrst-order analysis 
such as the present. In the work which follows, as in the 
calculations of Guderley and Toshihara, the supersonic 
region at the leading edge will therefore be disregarded. 

FORMULATION OF BOUNDARY-VALUE PROBLEM 

As is reference 2, the analysis is based on the equations of 
the transonic small-disturbance theory with the stream func- 
tion as the dependent variable. If the effects of the flow 
at the nose are ignored, the problem of the wedge at angle 
of attack a is then readily formulated as a boundary-value 
problem in the hodograph plane. To solve this problem for 
vanishingly small a, it will be assumed that the solution ip 
at angle of attack can be expressed as the sum of the basic 
solution Ip previously obtained at zero angle plus a perturba- 
tion term cop', where is a function which does not itself 
involve a. By consideration of the difference between the 
boundary-value problems for ip and ip, a problem for the 
perturbation fimction ip' can be formulated. “IJhe boundaries 
for this problem turn out to be the same as those for the 
problem at zero angle, and the boundary values themselves 
appear in terms of ip. The details of these matters will now 
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be giveu. The reader who is interested only in the results 
can proceed directly to the later section on Chordwise 
Distribution of Lift. 

Basic equations. — ^The basic equations will be taken in 
the form given in reference 2, that is, in terms of small dis- 
turbances from the critical speed The independent 

variables are the normalized speed ^ and the norm al ized 
inclination d as defined by the relations 


"’-VJa^-1 

(la) 

d-f 2 V ^ 

(lb) 


where 

y local speed of flow 

0 local inclination of flow relative to direction of free 
stream 

free-stream speed 

a* critical speed (i. e., speed at which the speed of flow 
and the speed of soimd are equal) 
y ratio of specific heats 

Use of these variables is equivalent to introducing the rules 
for transonic similarity. In terms of the foregoing hodograph 
variables, the differential equation for the stream fimction 
^ as given by the transonic small-disturbance theory is 

( 2 ) 

This is essentially the linear differential equation first studied 
by Tricomi (ref. 11). It is elliptic for ^<0 (subsonic speeds) 
and hyperbolic for 5>0 (supersonic speeds). 

The transformation from the hodograph to the physical 
plane is governed by the differential relations 

^ 1_ (3a) 

2 J 


tities rj and 6 are themselves the normalized quantities 
defined by equations (1). 

Problem at zero angle of attack. — ^When the angle of 
attack is zero, the localized region of supersonic flow at the 
leading edge disappears from figure 1, and the flow field 
becomes symmetrical about the chord line. The correspond- 
ing boundary-value problem in the plane has been set 
forth in reference 2. It is restated in figure 3, where both 
the upper and lower halves of the flow field are now included. 
In this representation, the surfaces of the wedge appear as 
the semi-infinite horizontal lines OB and OL, and the sub- 
sonic portion of the shock polar appears as the curve NAE. 


B 


const/ T}<irg{7],9„)cf7} 
^ J~a> 

B 

ijf ( 17 , ) » 0 

. 



Figubb 3. — Boundary-value problem at zero angle of attack in 

plane. 






where x=x(ij,d) and y=y{vfi) physical coordinates 
(horizontal and vertical, respectively), corresponding to a 
given velocity Tj,d. The symbol p+ denotes the fluid density at 
the critical speed a^. Within the approximation of the 
transonic small-disturbance theory, the pressure coefficient 
can be calculated from the relation 


If the stream function for zero angle of attack is denoted 
by ^=^'( 77 , 0 ), the differential equation to be satisfied here is 
given by equation (2) as 


( 6 ) 

The requirement that the flow shall be tangent to the surfaces 
of the wedge provides the boundary conditions 

i?(i7,±fl»)=0 for 77^0 (7) 


Gp=-2'^^^=-2(yja^-lKv-l) (4) 


The local Mach number is related to the speed of flow by the 
equation 


ikP-1 V 
y+1 a* 


(5) 


Eor simplicity of notation, the tilde will be omitted from 
the symbols ^ and 0 in the remainder of the development. 
It is to be understood, unless stated otherwise, that the quan- 


where d„ denotes the normalized half-angle at the leading 
edge. The stagnation point at the leading edge is represented 
in the present theory by the condition that 

^^0 for 77— CO, — Ba^O^Ou ( 8 ) 

Along the shock polar NAE, the relations for an oblique 
shock wave require that 

(9) 


* The theory could equally well be formulated In terms of disturbances from the free-stream speed Va (see, e. g., refo, 9 and 10). For a d i scussion of the relationship between this latter 
formolation and that used in the present paper, see reference 10. 
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for 

Along the sonic line, boundary conditions axe prescribed 
which represent the influence exerted on the subsonic field 
by the interdependent portion of the supersonic expansion 
fans. On the basis of the procedures given in reference 2, 
this influence can be represented completely b'y the require- 
ment that 


7,(0, «)+&£_ (10) 

where the upper signs apply for 1 ^0 and the lower signs 
for — The constant which appears here is 

given by 




Si/aps (1/3)' 


iO.3429 


where r(l/3) is. the gamma function of the axgument 1/3. 
The use of the relations (10) as boundary conditions along 
the sonic line reduces the transonic problem of the flow over 
the front wedge to a pm-ely eUip tic" problem in the hodograph 
plane. 

In addition to the foregoing conditions, a further condition 
is necessary to assure that the elution for ^ will give the 
proper scale when transformed to the physical plane. This 
is furnished, for example, by the following expression for the 
half-chord of the profile, found by integrating equation (3a) 
over either OB or OL: 


2 

2 


j~ (T+l) (VJ(h~ 


• 1 ) 


- 1/3 J*0 


V'f'e (v,±Sa) dn 


( 11 ) 


If the chord of the profile is given, this condition, together 
with the previous conditions (7) through (10),issuflScientto 
determine a uilique solution to the problem. 

It is obvious from the nature of the boimdary-value prob- 
lem (and also from considerations of symmetry in- the phy- 
sical plane) that the solution for must be antisymmetric 
with respect to 6. The problem can be simplified, there- 
fore, by discarding the lower half of the hodograph and 
replacing it by the condition 


^(17, 0)=0 for 17^-1 (12) 

The resulting problem is readily solved with numerical 
methods by assuming an arbitrary value of ^ at some point 
(as, e. g., point E), sol-dng for ^ in the upper half of the 
hodograph subject to the conditions (7), (8), (9), (10), and 
(12), and then adjusting the solution to satisfy condition (11). 

It is apparent from the boundary conditions that the 
solution of the foregoing problem -will depend on the value 
of the parameter which defines the position of the upper 
and lower boimdories in the hodograph. This parameter is 


related to the more familiar transonic similarity parameter 
by the relation 

^ [(T-l-l) 

where tjc is the thickness ratio of the complete double-wedge 
profile.® In reference 2 the solution of the foregoing problem 
has been carried out for four values of 

Problem at angle of attack. — ^If the supersonic region at the 
leading edge is ignored, the boundary-value problem for 
the wedge at angle of attack appears in the rj,d plane as 
sho-wn in figure 4. The primary difference between this and 
the pre-rious figure is that the lines OB and OL, which repre- 
sent the surfaces of the wedge, have each been displaced 
downward by an amount a, where a is the angle of attack 
normalized in the same manner as the other angles of inclina- 
tion (cf. eq. (lb) ) .* 


e 



Figuhb 4. — Boundary-value problem at angle of attack in jj,ff plane. 


The stream function at angle of attack -will be denoted 
here by ^=^(i7,0;«), the latter notation being used to indi- 
cate the dependence of ^ upon the parameter a. The 
function ^ must satisfy the differential equation (2), which 
is now -written. 

'f'in~^vfee=0 (14) 

The boundary conditions at the surface of the wedge now 
require that 

^(»7,±0ui — o;;a)=0 for 17 ^0 (15) 

while the condition at the leading edge becomes 

^—>0 for 57 -h> — ooj — da — a^B^Ba — a (16) 


> Tbrougbout the present report, the expressUms Ibr the transonic similarity variables have been altered from those used In the earlier accotmt of this wort (see footnote 1) by replacing the 
qaanHty( 7 +l)v'hereverltappearslnthesevarIabIesby the product ( 7 +l)ArcD’- This procednre Is suggested by the derivation of the similarity rules given by Bnse m an n In reference 9. It 
has been shown by Sprelter (refs. 10 and 12) In several spedfro examples to Improve the accuracy of the smaB'dlsturbanoe theory as compared both with experiment and with exact theory. 

r In reference 8, Quderley and Yoshlhara find It convenient to obtain the angle of attack by" holding the profile fixed and dianglng the Inclination of the free stream. This procedure, k 
applied In the present case, would reqalre the eventual calculatlou of the second derivatlvosof^ ou the shock polar. The present procedure, which holds the free stream fixed and changes the- 
attitudeoftheprofile, requires the calculation of only a first derivative of f at the surface of the wedge. Since the accuracy of numerical differentiation decreases with Increasing order of the- 
derivative, the present approach Is to be preferred In a numedcal analysis. 
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The shock polar NAE is imaltered from the previous prob- 
lem, and the condition on this boundary has the same form 
as before. The conditions along the segments BE and LN 
of the sonic line are now 




(17) 


± Ou— a) =Hv, ± 0b) — ccil/eiv, ± 0») + 0 (a*) (22a) 

(i7)±0o) — 9(’?,±0 ui)+O(q:^ (22b) 

B these expansions are substituted into equation (21) and 
'f'ivjdzOa) set equal to zero in accord with the boimdary con- 
dition (7), one obtains finally for vanishingly small a 


where the lower limit of the integral has been changed in 
accord with the displacement of the points B and L. The 
upper signs in equation (17) now apply for l^d^da—a 
and the lower signs for — — 1. An expression 
for the half-chord of the profile can be found again byinte- 
gratmg equation (3a) over the line OB or OL, which givte 

^^[ev±MMJ-£ as, 

If the chord of the airfoil is specified — say the same as at 
zero angle of attack — then the foregoing conditions are 
sufficient to determine a solution. No simplification based 
on symmetry considerations is possible in the present case. 

Perturbation problem. — The problem of the preceding 
section could conceivably be solved by numerical methods — 
though with great labor — ^for arbitrary values of a. Efforts 
in this direction would hardly be justified, however, in*view 
of the fundamental omission of the localized supersonic 
flow at the leading edge. It is more reasonable to examine 
the problem for vanishingly small a, where this ommission 
is valid and where' there is hope that the amount of labor 
might be reduced. 

To proceed along these lines, it is assumed that yf/(T],d;a) 
may be expanded in a power series of the form 

HVj^',<^)=Hv,S'fi)+ayf'a(v,S',0)+0 (a*) 

where, for present purposes, only terms to order a need be 
retained. The first term on the right represents the solution 
at a=0 and is thus identical with the function ^(17,0) pre- 
viously introduced. The second term will be abbreviated 
by means of the notation If terms of 

0(a*) are discarded, the expression for ^ can then be written 

Hv,6]cc)=^{i],0)+a ^'(v,0) (19) 

By comparison of the previous boimdary-value problems 
for f and a problem for the perturbation function }{/,- 
will now be formulated. 

The differential equation for ^ follows at once from 
the differential equations (6) and (14) and the substitution 
(19). It is obviously of the same form as the previous 
equations, that is, 

^p'es=0 (20) 


r(.v,±e.)=Mv,±e.) for (23) 

This is the boundary condition for appropriate to the 
surface of the profile. It will be noted that the condition is 
applied in the hodograph at the original, imdisploced loca- 
tion of the surface (i. e., 6= ± 0„) . The condition depends for 
its apphcation on a knowledge of the basic solution \f>. 

The boundary condition for yp' at the leading edge follows 
directly from the conditions (8) and (16). It is the same as 
the corresponding condition for p, that'is, 

^'^0 for i 7 ->— 00 ,— (24) 

As was indicated, the functions ^ and p both satisfy the 
same linear; homogeneous boundary condition on the shock 
polar. It follows, as in the case of the differential equation, 
that the condition for p' on the polar is again the same, that 
is, 

(26) 

for 

0=±(i — >?) Vi+’/j 

The treatment of the boimdaiy condition along the sonic 
fine is complicated by the fact that the parameter a appears 
in the condition (17) as a term in the lower limit of the inte- 
gral. For simplicity, the details will be confined here to the 
upper segment BE of the sonic liue. For this segment, 
condition (17) becomes, after substitution from equation (19), 


^»(O>0)+’a^ 7(0,0)+^ 






applicable for lS6-^d„~a. To simplify this equation, the 
first integral is rewritten 


r MOA) 

j9.-a(0l-0)*" 


d6i= 


^ UOA) 

i.(0i-0)*/® 


d0i— 


^«(O)0i) 


dBi 


(27) 


It can be shown from Guderley’s analysis of flow at a convex 
comer (ref. 13) that, for vanishingly small values of (0»— 0), 
the variation of ip along the, sonic fine must be of the form 


The boundary conditions appropriate to the surface of 
the wedge are established as follows: The boundary condi- 
tion (15) for }p is first rewritten, with the aid of the substitu- 
tion (19), in the form 


rp(o,e)=0(.e^-e)*/^ (28) 

where G m & constant for any given value of 0„. Differen- 
tiating this relation, one obtains 


j{v,±e„-^cc) + ccriv,±d„~a)=0 ( 21 ) 


?9(O,0i)~(0„-0i)>« 


By expanding in Taylor’s series about the fines 0=±0„, the 
functions p and p' can be written 


Substitution of this result into the second integral' of equa- 
tion (27) yields the fact that this integral must be proper- 
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tional to a*^. The first integral in equation (26) can thus 
be Avritten 


MOA) 





dd,-\-Q{a.^^ 


(29) 


The second integral of equation (26) can be treated shnilarly 
by first rewriting it as 





tA^(0,gi) 



%(0,ei). 


(30) 


To deduce the variation of rp' for vanishingly small (6a— d), 
it is first noted that a result s imil ar to equation (28) must 
also hold for the variation of \}/ relative to the displaced 
location of the shoulder, that is, 

rpi0,d)=C(6a-a-d)^ 


The quantity C=C(a) is a differentiable function of a 
which reduces to C when a= 0 . Since a will eventually be 
made less than any assignable value of (d„—6), this expres- 
sion' may be expanded in the form 

^O,fl)=(7(0«-0)^'^[l-|^+O(a*)] (31) 


It remains to impose the condition that the chord of the 
airfoil must remain rmaltered during change in angle of 
attack. To express this condition in terms of equation 
(19) is first substituted into equation (18) to obtain 

2 r (-Y-H)(7„K-l) T/^ 

2 L 2 J 

r I7[^9(i7,±0«— a)+o:i/^9(i7,db0u,— a)]dij (35) 

C/ — CO 

As in the treatment of the boundary conditions along the 
upper and lower boundaries, Taylor’s expansion gives 

^«(i?,±0» — oc)=i/e(v, ±0tr ) — ccpeeivi i^irl+Ofa®) (36a) 

'l''e(‘n, ±0tt)— a^'«(i?, ±0 »)+O(q:*) (36b) 


It can be inferred dhectly from the boundary-value problem 
for \p that ^ 99 ( 1 ?, ±0»)=O, so that the term involving this 
quantity may be dropped from equation (36a) . Substitution 
of equations (36) into equation (35) and application of the 
previous expression ( 11 ) leads, for vanishingly small a, to the 
condition that 


j_^V'l''s(v,±Si^dv=0 (37) 


Now it follows from the definition of yf/ that 

«-*0 ^ 

Substitution from equations (28) and (31) thus gives for 
the variation of \p' in the vicinity of the shoulder 

no, 6) =lim r ^ (0„-e)«’-| O(0„-0)'/’-bO(a)1 
a -»0 L ® ^ J 

or 

no,B)=C'(da-e)^^-^ (7(0„-0)'« (32) 


where O'^Ca(0). This means that for vanishingly small 

(d,.-0) 

' 4''e(o,ed~(e.-6i)-^'^ ■ 


On the basis of this result equation (30) can be written 



^'ei0,dy) 

(0i-0)*/3 



(0i-0)^« 


ddi+0(a^'^ 


(33) 


If equations (29) and (33) are substituted into equation (26) 
and the boimdary condition ( 10 ) is taken into account, one 
then obtains for vanishingly small a 




r» mo, 61 ) 
Jo, (6:-dr^ 


ddi=0 


(34) 


where l^B^da. The boundary condition for yp' along the 
upper segment of the sonic line is thus the same as the con- 
dition for yp. The same result can be shown to hold along 
the lower segment. 


The boimdaiy conditions (23), (24), (25), (34), and (37) are 
sufificient to determine the solution for yp' in the hodograph. 

As with yp, the boimdary-value problem for yp' can be 
simplified from considerations of symmetry. Since yp is 
antisymmetric •with respect to 6, the nonhomogeneous 
boundary condition (23) which is imposed on yp' along the 
upper and lower boundaries must be symmetric in thi.q 
variable. The remaining conditions, which are all homo- 
geneous, are also synometric. It follows that yp' itself must 
be a symmetric function of 0 .® The problem can therefore 
be simplified by again eliminating the lower half of the 
hodograph and substituting in this case the condition that 

'A's(’?,'0)=0for 1 (38) 

The problem which is finally to be solved is thus as sum- 
marized in figure 5. The boundaries for this problem are 


B 



• This result eon also be argued directly from considerations In the physical plane. It is necessary to mate two observations as lOHows: (1) Since the profile Itsell Is symmetric about the 
chord Hue, the flow field at a negative angle of attach must bo the Inverted Image 0 / the flow field at an equal poslUve angle. (2) To be consistent with the basic perturbation assumption. It 
must bo presumed that all changes In the flow field are smooth funcUons of angle of attach at a— 0. These statements taten together Imply that the vertical distance between any two points 
of equal 1 ) and corresponding positive and negative 9 Is, to a first order, unaftectod by angle of sttaok. It follows that, for sufllolently small ct, the Increments in 5 and ^ between the two polnta 
are equal and hence, on the basis of equation (19), that the value of at the two points la the same. 
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identical with those used to obtain yj/. The boundary con- 
ditions are also identical insofar as the shock polar and sonic 
line are concerned. The only differences between the two 
problems are in the conditions imposed along the boundaries 
OB and OA. As was the case with tj/, the solution here must 
be a fimction of 

Because of the nature of the integral condition along the 
upper boimdary OB, a direct solution for yp' is not feasible by 
numerical methods. To obtain a solution, therefore, the 
problem is broken down into two subsidiary problems by 
means of the substitution 

P'=P'A.+ hyp'B (39) 

where 6 is a constant whose value is to be determined. 
Boundary-value problems for yp'j. and yp's are-then defined 
as shown in figme 6. In both problems the integral condi- 




Figubb 6. — Subsidiary perturbation problems. 

tion along the upper boundary OB is ignored for the time 
being, and an arbitrary condition is introduced instead at 
the point E. In the problem for yp'j,, only the nonhomoge- 
neous condition (23) is imposed at the upper boimdary, and 
the condition at E is the homogeneous one that ^a= 0. 
In the problem for yj/'a, the homogeneous condition ^'b= 0 is 
imposed along OB, and the condition at E is that yp's has an 
arbitrary value The conditions at the remaining 

boundaries are the same as in figure 5 and are therefore not 
repeated here. It is apparent that a superposition of 
and yp'a will constitute a solution of the original problem pro- 
vided the value of 6 is adjusted so that the integral condition 
(37) is satisfied on the upper boundary. The necessary 
equation for h is found by substituting the expression (39) 
into condition (37) and is 


J 

[* WpAfiViBw) dt] 

] 

f vpBi(.v,Ba) d-n 

“• eo 


Eolations for quantities in physical plane. — ^To complete 
the fundamental analysis, relations must be established 
between yp' and the relevant quantities in the physical plane. 
Let x=x(i], 0) and j7=^(Vj®) denote the coordinates at which 
a given velocity rj,d is found m the physical plane when the 
profile is at zero an^e of attack. As shown in reference 2 
(pp. 22-23), the transformation equations (3), when apphed 
to the case of zero angle of attack (and written in the present 
notation), can be put in the generalized dimensionless form ® 


d(^ 

41a 


(Ala) 

[(y+i)MimV'^d(^= 

(p,dr,+p,de)= 

Ala 

Ala 





(Alb) 

where I„ represents the integral 





j: 


V^e(v,0v) dv 


(42) 


By taking the origin of the physical coordi^tes at the lead- 
ing edge and introducing the notation X=$lo and Y= 
[(rd-l)Af*(t/c)]'/®(^/c), equations (41) can be integated to 
give 


^(‘n,0)=‘-^Jj2vPedv+p^de) (43al 


nv,e) 


(26g^- 

4/„ 


(43b) 


The integration in equation (43a) is performed in the hodo- 
graph over any curve C which begins at rj=— oo and ends at 
the point r],d. The generalized coordinates X and Y at 
which the same velocity rj,6 is to be found when the airfoil 
is at angle of attack are given correspondingly by 


X(ij,e;a)=^ Jj2v'f'idv+'f'^) 

W n. ^ (2^.r)’'® , 

Y{it,e-,a)=^^-^^- yp 


(44a) 

(44b) 


The integration in equation (44a) is considered to be taken 
over the same cmve C as before.'^ The integral is now 
given by 



(46) 


It can be shown from equations (19), (36), and (37) that for 
vanishingly small a 

!»=/« (46) 

. Equations (44) can now be specialized in the light of the 
basic perturbation assumption. This assumption imphes 
at once that the coordinates X and Y in the physical plane 
must be expressible in the form 

X(jjfi-a)=X{-n,e)+aX' (47a) 

Y{-n,d-,a)=Y{-n,B)-^aY'{ii,B) (47b) 

where X'{ri,e)=Xair},B-fi) and Y'{ri,B)^Ya{ri>B-,0). If ex- 
pressions (19) and (47) are substituted into equations (44), 

and equations (43) and (46) are taken into account, the 
following relations are finally obtained for X' and Y' in 
terms of p': 

r {2i]p'edi]+p^ (48a) 

41a J a 

yp' (48b) 


• See footnote 3. 

I If C lies slightly ontsWo the domain In whldi Is defined— as will bo the case, for example, when the Integration Is taken over the upper surface of the wedge In Its undisplaced position— 
Is to be thought of as being oontlnned analytically outside the bonndary. 
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The foregoing equations (48) give the initial rate of move- 
ment with angle of attack of a point of fixed velocity i],0- . 
One requires for practical application, however, the rate of 
change of tj and 0 at a point of fixed location X, Y. Equations 
relating the two sets of derivatives can he obtained as 
follows: If and 6 are regarded in the physical plane as 
functions of X,Y, and a— that is, i]—ri(X,Y;a) and 
0 = 0 (X,F;o:)— then the corresponding total differentials are 

dr]—rixdX-\-r]YdY-\-rjada (49a) 

dd=6xdX-{-9YdY-\-8ada (49b) 

Consistent with the basic perturbation assumption, and 8 

can bo written 

v{X,Y;a)=v{X,Y}+c^'(X,Y) (50a) 

d{X,Y-,a)=e{X,Y)+aB'{X,Y) (50b) 

where Ij and 6 represent the conditions at a given point X, Y 

at zero angle of attack _md 17 ' and d' are defined by 
r,'(X,T)mr,,(X,Y;0), 6'(X,Y)=d.(X,Y;0). In view of equa- 
tions (60), equations (49) can be written for vanishingly 
small a 

di]='ijxdX-^rjydY-\-i]'da (51a) 

d6=^6^X-\-0ydY-\-6'dci (51b) 

Similarly, from equations (47), one can write for the differ- 

entials of X and Y as functions of 17 , 6, and a 

dX=X4ri+Xede+X’da 

dY=Y^^+Yedd+Y'da 

from which 

X4v+Xede=dX-X'da (52a) 

Y^t]-\-Yedd=dY—Y'da (52b) 


graph. The coordinates at which these derivatives apply 
axe found from the solution at zero angle of attack by means 
of equations (43). 

The foregoing equations are considerably simplified when 
applied at the surface of a wedge profile. Here the boundary 
condition is that ^ is constant on a line of constant 9 (cf. eq. 
(7)), with the result that Equation (54a), for ex- 

ample, can thus be written as simply 


,'(X,±0) (S6) 

where the upper signs pertain to the upper surface and the 
lower signs to the lower surface. The corresponding rate of 
change of pressure coefficient is found by differentiating 
equation (4) with respect to angle of attack. If a is used 
now to denote the true angle of attack (related to the previ- 
ously used, normalized angle of attack by an equation like 
(lb)), such differentiation then gives 



-2 


' 2 1 

J+ 1 / 


V 


/ 


Here 17 ' is stiff the derivative with respect to the normalized 
angle as given by equations (54a) or (55). With the aid of 
equations (5) and (13), this result can be rewritten 


[{y+l)MJmr^ 2 ( 28 „) i / 3 , 7 ' ( 66 ) 


It can be seen from equation (56) and the symmetry prop- 
erties of ^ and ^ that ij' must be of equal magnitude but 
opposite sign on the upper and lower surfaces of the profile. 
If the local lift coefficient is represented by = (piower — 

2^npper)/2„) it then follows from equation (56) that 

=4(28„)*''’7,'(X-1-0) (57) 


Solution of equations (62) for dr] and dd and comparison of 
the results with the alternative expressions (51) gives finally 
for 17 ' and 9' 


where the notation 17 ' (Z", - 1 - 0 ) indicates that the value is to 
be taken on the upper surface of the profile. Substitution 
from equation (65) gives finally 


, 

Xje-XeT, 

T^'-X,Y' 

X,7e-TeT, 


(53a) 

(53b) 


These equations can be put in more fcectly useful form by 
evaluating the derivatives of H and Y from equations (43) 
and substituting for X' and Y' from equations (48). There 
results finally 


l(y+i)MJit/c)y>’‘ 


r d{^vlg.J^ _ 

4(28„)>« r 

da Jc=o 

vl'eivjSa) J- 


V'l''e(j],9w)dv 

(58) 


By means of this equation, the initial rate of growth of lift 
at any chordwise station can be obtained. Since ^ and 
axe both fimctions of the parameter 9a, the generalized 
quantity which appears on the left-hand side of equation 
(67) is also a function of this parameter. These results are 
in conformity with the rules for transonic similarity (see, 
e. g., ref. 10 ). 


[i}ejj.2vf'edv+r,d9)-^,r] (54a) 

9'(XyY)— ^^^ - _^i j^?',J^(2i?'/^edi7-l-iA'i^)— (54b) 

By means of these equations the initial rate of change of 
7} and 9 at some fixed point in the physical plane can be' 
calculated corresponding to any chosen location in the hodo- 


MBTHOD OF SOLUTION 

As in the previous calculations of the boundary-value 
problems for and iI^'b can be solved through the use of 
finite-difference equations and relaxation techniques. A 
detailed description of the general method has been given in 
reference 2 and need not be repeated here. Most of the 
necessary finite-difference equations — ^notably the tedious 
ones along the shock polar and sonic line — can be taken over 
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e. g.jfig. 27 of ref. 2). Staxting from the kno\vn solution for 
one first employs equations (48) to compute the initial 
rate of movement of the points at which the Mach lines of 
the basic characteristics net meet the sonic line. Using these 
results and the known slope of the segments of the basic not, 
one then proceeds stepwise along consecutive downgoing 
Mach lines, calculating the initial rate of movement of suc- 
cessive intersection points on each line. By application of 
the proper boimdary conditions at the surface of the wedge 
the value of at the surface is finally determined, and from 
this the initial distribution of lift is calculated. The details 
of the procedure are given in Appendix B. 

RESULTS AND DISCUSSION 

Calculations of the lift have been carried out, foUoAving 
the methods just outlined, for the same values of used in 
the work at zero lift, namely, 1.3, 1.6, 2.4, and 4.2. These 
valnes correspond, respectively (see eq. (13)), to values of 
the similarity parameter of 1.058, 0.921, 0.703, and 0.484. 

To illustrate the results for the front wedge in the hodo- 
graph, figures 7 to 9 have been prepared showing the vari- 
ation of yf/'j,, rp's, and for 0„=1.6®. The results for 


* For the calcolatloiiof^'B In. tlUs example, use was made of 236 lattice points dlstrlbated as shown for ^^Infigore 22 of re&rence 2. For ^.i,3S0 points were used with a distribution appro* 
priato to the altered behavior of the dependent variable. 



directly from the previous work. The only equations which 
need be altered are those directly infiuenced by the change in 
boundary conditions on the upper bmmdary and on the 
horizontal axis. The only real difficulty from this source is 
encountered in the solntion for in the vicinity of the 
shoulder (point B in fig. 6). At the shoulder itself, the 
boundary conditions require a singularity in the first deriva- 
tives of which means that any purely numerical treat- 
ment would be of doubtful vahdity in the vicinity of this 
point. This difficulty is overcome by subtracting out an 
analytical solution of the proper singular form and then 
working locally with the difiference between this solution 
and the desired unknown. The singular solution is obtained 
from the general results of Guderley (ref. 13) and is expressed 
in terms of hypergeometric functions. The details of this 
and other matters regarding the numerical calculations for 
the front half of the profile are given in Appendix A. 

With the solution known for the front half of the profile, 
the calculation of the lift on the rear half is a simple matter. 
The computations are carried out in the physical plane and 
are based on the characteristics net previously constructed 
for the flow over the rear wedge at zero angle of attack (see. 
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'i''a (fig- 7) jure only slightly different from those previously 
shoivn for tp in figure 23 of reference 2. As before, a rapid, 
(but regular) variation is apparent in the dependent varia- 
ble in the vicinity of the point 17=0, 6=1. The results for 
p'j, (fig. 8) show a rapid variation near the point i;=0, 
6= 6 a. This is a consequence of the previously mentioned 
singularity in the first derivatives of at that point. The 
values of p' (fig. 9) are found in the present case from the 
equation 0.5348 \P'b (cf. eq. (39)). They exhibit 

the same behavior as does in the vicinity of the singular 
point but differ markedly in other parts of the field. For 
reference, the numerical values from which figures 7 and 8 
were plotted are given m tabular form at the end of the 
report. 

The complete results for the lift of the profile are given in 
figures 10 through 13. These results will be discussed in the 
following paragraphs. 

CHOBDWISE DISTHIBOnON OF LIFT 

Figme 10 shows the calculated lift distribution, in transonic 
similarity form, for the four values of considered in the 
present work. Also included are the results for f„=0 
(M„=l) given by Guderley and Yoshihara in reference 8. 
As can be seen, the same general features are apparent for 
aU values of the similarity parameter. 

At the leading edge of the profile, for example, the calcu- 
lated lift tends in all cases toward infinity. This physically 
impossible result is, of course, well known from the linear 
theory of airfoils at subsonic speeds. It is a result of the 
obvious failure of the small-disturbance approximations to 
conform with the actual phenomena in the vicinity of the 
leading edge. This local fadme of the theory is faiown in 
the linear, subsonic case to be of httle consequence insofar 


as the over-all lift is concerned. It may be presumed that 
a s imil ar situation exists here. 

As one proceeds rearward from the leading edge, the lift 
distribution falls more or less rapidly, reaching a value of zero 
directly forward of the shoulder. This latter result could have 
been foreseen, since the speed on both the upper and lower 
surfaces has a fixed (i. e., sonic) value at this location. 
Directly to the rear of the shoulder, the lift distribution 
starts anew from zero. This must obviously he the case, 
since the expansion from sonic speed is, in Prandtl-Meyer flow, 
a unique fimction of the local turning angle, which is the same 
for both surfaces. Rearward from the shoiflder the lift 
increases monotonically to a relatively small, finite value 
at the trailing edge. 

OvOT the front wedge, the four ciuves of the present study 
exhibit a uniform progression with respect to The curve 
of Guderley and Yoshihara, however, crosses the present 
curves at several points. The reasons for this are not clear, 
though it is hi^y unlik ely that such a result eould be in fact 
correct. The observed behavior may be due to sbme con- 
sistent inaccuracy in the present numerical approach or to 
the approximations introduced by Guderley and Yoshfliara 
in satisfying the boundary conditions for the interdependent 
portion of the supersonic expansion fan. Over the rear wedge, 
the present computations give virtually'a single curve for the 
four values of There is again, however, a small incon- 
sistency with the results given by Guderley and Yoshihara. 
This is as might be expected if the calculated flow over the 
front wedge is in error in either case. 

UFT-CUBVE SLOPE 

Figme 11 shows the generalized slope of the lift cmwo at 
zero angle of attack plotted as a function of the transonic 


O 

a 
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Fiqtob 11. — Generalized lift-curve slope as a function of transonic similarity parameter. 


similarity parameter. Results obtained on the basis of the 
transonic small-disturbance theory are shown by three 
solid-line curves. Each of these curves consists of two 
segments separated by a gap within which the curve cannot 
be defined on the basis of the available results. The upper- 
most of the three curves gives the lift of the complete profile; 
the other two show the division of lift between the front 
and rear wedges. 

The left-hand segment of each of the cmwes in figure 11 
shows the variation of Uft-curve slope over most of the 
range of flight speed in which the bow wave is de.tached, 
which is the range of primary concern in the present analy- 
sis. The calculated points from which these curves were 
di'awn are shown in the figure. The points denoted by 
squares were obtained by mechanical integration of the lift- 
distribution curves of figure 10.® The circled points on the 
vertical axis were located on the basis of the work of Guderley 
and Toshihara.*® 

The right-hand segment of the curves in figme 11 shows 
the variation of lift-curve slope in the range of flight speed 
in which the bow wave is attached and the flow is com- 
pletely supersonic. To the order of accuracy of the present 
theory, this condition exists for the double-wedge profile at 
zero angle of attack when ^2‘^=1.260.“ Above this 
value, results completely consistent with the fundamental 


assumptions of the transonic small-disturbance theory can 
easily be obtained by analytical methods. To this end, one 
need only presume that the speed is constant on each straight- 
line portion of the airfoil surface, a condition which is actually 
fulfilled over most of the pertinent range of The neces- 
sary procedures are outlined in Appendix C. To the accu- 
racy of the transonic small-disturbance theory, the results 
provide an exact solution for the lift-curve slope of the front 
wedge for all values of in the range of completely super- 
sonic flow. For the rear wedge — and hence for the com- 
plete profile — the solution is exact down to a limiting value 
of somewhat greater than 1.260. Below this limit the 
intei’action of the shock wave from the bow and the expan- 
sion fan from the shoulder influences the flow over the rear 
wedge, with the result that the condition of constant speed 
is not satisfied. The position of this limit is difficult to 
determine exactly. As shown in Appendix C, however, it 
must he at a value of less than 1.287. The curves for 
the rear wedge and complete profile are thus approximate 
for at least a portion of the interval from 1.287 to 1.260 and 
are therefore shown dotted in this range. It can be demon- 
strated that inclusion of the interaction effects in the analy- 
sis woifld cause an increase in the computed lift for the rear 
wedge. Exact results would thus he somewhere above the 
dotted portion of the curves in figure 11. 


• As In the earlier calculntlonsof the drjgcoefflclontat zero angle (cf.pp.9and24of rof.^, the Integration over a smaU Interval near the leading edge was carried out analytically on the basis 
of an asymptoUo ropresSntatlon of tho solution In the hodograph plane. 

™ In flgnre 11, as In figure 13 that follows, the ordinates for fa>“ LOSS have been changed somewhat from the valnos given In an earlier publication (soe footnote 1). This was done to correct 
an error made In the original Integration of the lift dlstributlou on tho front wedge. 

>• Attachment of the wave takes place at the somewhat lower value of $a> =3/(4)V«— L19L For 1.191<f co <1.280 the wave Is attached but the flow behind It Is still subsonic. 
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The most interesting aspect of figure 11 is the behavior 
of the lift in the vicinity of shock attachment. Despite 
the gap in the curves in this vicinity, it is obvious that the 
lift-curve slope of the complete profile must attain a maxi- 
mum somewhere in the range from {„= 1.058 to 1.287. 
This is in marked — and somewhat surprising — contrast to 
the previous- results-fop the drag coefficient at zero angle of 
attack, which was found (ref. 2) to decrease monotonically 
as the similarity parameter increased above zero. The peak 
in the curve in the present case is accompanied by a similar 
variation in the lift-curve slope of the front wedge. The 
results for the rear wedge may or may not pass through a 
minimum in the same range of 

A determination of the exact shape of the curves in the 
vicinity of shock attachment is not feasible on the basis of 
the present laborious methods. The existing curve for the 
complete profile does show a maximum in the range of com- 
pletely supersonic flow, but this is in the portion of the 
range in which the computed curve is known to be errone- 
ously low. If exact results were available for all values of 
the maximum would undoubtedly be somewhat higher 
and displaced somewhat to the left. The infinity which 
appears in the slope of the curve at {„=1 .260 (see Appendix 
C) would probably disappear as well. The lift of the rear 
wedge, which now ^es to zero at f„=1.260, would pre- 
sumably remain finite throughout. 

Tflthin the transonic range itself, the curves of figure 11 
show little variation for some distance above a similarity 
parameter of zero. This is in accord with Guderley^s ma- 
ly tical study of two-dimensional flows with a free-stream 
Mach number dose to 1 (ref. 14). Guderley’s work shows, 
in particular, that the curves of flgure 11 should have zero 
slope at f„=0. The figure has been drawn to conform 
with this requirement. It appears from the present work 
that, for practical purposes, Guderley’s result may be taken 
as valid in the range of up to about 0.5. The same 
result was found in reference 2 with regard to the drag 
coeffident of the complete profile at zero angle of attack. 

Over most of the range of completdy supersonic flow, the 
lift-curve slope of the complete profile exhibits the type of 
variation well known from linear theory. This latter theory 
gives for the lift-curve slope of all thin profiles 


dci 4 


(59) 


which can be written in terms of the transonic similarity 
variables as 


[(T-fl)M„*(t/c)]>«g=|^ (60) 


The dashed curve in figure 11 is based on this equation. 
There is' considerable quantitative difference between the 
linear and nonlinear results for values of just above 
1.287. As increases, however, the curves given by the - 
two theories appear to converge. This latter behavior is in 
accord with Spreiter’s considerations (ref. 10) regarding the 
basic relationship between the linear and nonlinear theories. 

To put the results in more familiar form, the lift-curve 
slope of the complete profile has been replotted in figure 12 
as a function of Mach number for 'y=1.4. The results of 



Free-stream Mach number, /Wq) 


Figure 12. — Lift-curve slope as a function of Mach number for sovoral 
thickness ratios (7=1.4) 


linear theory give a unique curve defined by equation (69). 
The no nlin ear, transonic theoiy provides a family of curves 
with thickness ratio as a parameter. As would be expected, 
the range of Mach nmnbers over which the linear thooiy is 
a poor approximation becomes smaller as the thickness ratio 
is reduced. It can be reasoned, in fact, that the nonlinear 
results must tend toward the results of the linear theory as 
f/c->0. 

CENTER OF UFT 

Figure 13 shows the chordwise position of the center of 
lift {x!c)i as a function of the transonic similarity parameter. 
The arrangement of the figme parallels that of figure 11, 
As before, the indicated points were calculated on the basis 
of the lift distributions of figure 10. The curve in the range 
of completely supersonic flow (f„^ 1.260) was obtained by 
means of the equations of Appendix 0. Only results for the 
complete profile are shown. 

The movement of the center of lift with increasing is of 
some interest. At f„=0(M„=l), the results of Guderley 
and Toshihara indicate a position about 29 percent of the 
chord aft of the leading edge. As the value of is increased, 

the center of lift first moves forward, slowly in the initial 
stages and then more i-apidly as the condition for shook 
attachment is approached. In the completely supeisonic 
range, this trend is reversed; the center of lift then moves 
aft toward the midchord location given by linear thooiy. 
Apparently, the reversal of the direction of motion must 
take place rather suddenly in the vicinity of shock attach- 
ment. The limit of forward movement cannot be specified, 
except to say that it must lie somewhere ahead of 25 percent 
of the chord (and probably aft of the leading edge). The 
dotted (i. e., inexact) portion of the curve passes precisely 
through the quarter-chord point at 1.260. (The cor- 
responding lift distribution is one of uniform lift on the front 
wedge and zero lift on the rear.) Because of the interaction 
effects previously discussed, an exact result would lie some- 
what above the dotted curve. 

CONCLUDING REMARKS . 

The present calculations add support to the growing con- 
clusion (see refs. 2, 5, 6, 7, and 14) that no marked changes 
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tako place in characteristics of airfoil sections as the free- 
stream Mach niunber passes through 1. The establishment 
of this conclusion must be regarded, in fact, as one of the 
major successes of recent research in transonic flow. In the 
present case, as in the previous study of the drag coefficient 
at zero lift, the variation of the aerodynamic quantities with 
freo-stream Mach number is most rapid in the vicinity of 
shock attachment. Unlike the behavior of the drag co- 
efficient, however, the variations here are large and charac- 
terized by a sudden reversal in the sign of the derivative. 
In drawing conclusions from these results it must be remem- 
bered, of course, that the theory assumes an inviscid medium 
and an airfoil of small thickness and infinite span. It also 
assumes, in effect, that at a given Mach munber the angle of 


attack is small compared with the difference between the 
actual wedge angle and the wedge angle that woifld provide 
shock attachment at zero lift. To what extent the results 
will be valid for viscous flows about finite-span airfoils at 
practically usable values of the thickness ratio and angle of 
attack is difficult to say. The effects of finite span, for 
example, will surely cause a reduction in the variations near 
shock attachment. In the present state of theoretical 
development, the study of these effects is a task for experi- 
ment. 

Ames Abeonautical Laboratohy 

National Advisory Committee por Abeonatjtics 
Moffett Field, Calif., Aug. 1, 1952 




APPENDIX A 

SOLUTION OP BOUNDARY-VALUE PROBLEM FOB FRONT WEDGE IN HODOGRAPH PLANE 


The solution of the boundary-value problems for and 
\}/'b was accomplished by finite-difference methods similar 
to those developed for the calculation of ^ in reference 2. The 
description here will be limited to the few featines wherein 
the present work departs from that discussed in the earlier 
paper. (See general remarks imder AIETHOD OF SOLU- 
TION.) The notation and sketches follow the conventions 
used in reference 2. 

FINITE-DIFFEEENCE EQUATIONS COMMON TO BOTH PROBLEMS 

The only finite-difference equations common to the prob- 
lems for and yj/s but not found in the problem for yf/ 
derive from the boimdary condition on the horizontal axis 
(see figs. 5 and 6). This condition is given for both problems 
by equation (38) and is iA'«(i7,0)=0 for 1. In the 
previous work, the finite-difference equations for lattice 
points located on a boundary were obtained by approxima- 
tion to the boundary condition itself. In the present case, 
the approximation to the differential equation will be 
employed, and the boundary condition incorporated through 
use of the equivalent symmetry property. 

Consider a typical point 0 on the horizontal axis as shown 
in figure 14. Point 3 is a fictitious lattice point located 

I 


A 


3 

Figure 14. — Point on horizontal axis. 

below the horizontal axis at 6=—A, where A is the lattice 
interval. The finite-difference approximation to the differ- 
ential equation (20) of the present text is given by equation 
(37) of reference 2 as 

r2+yp\-2vo{yf'\+rz) -2 (1 -2,o) ro=0 (Al) 

where rjo is the absicissa of point 0. The symmetry property- 
leading to the boundary condition (38) requires that 'P'z='P'i, 
so that for points on the horizontal axis equation (Al) 
reduces to 

\(''2+iA'4-Wi-2(l-2i7o)i^o=0 (A2) 

The point at the intersection of the horizontal axis and the 
shock polar needs special consideration. Figure 15 shows 
562 
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Figure 15. — Point at intersection of horizontal axis and shook polar. 

conditions at this point. Here, as before, point 4 is a 
fictitious point located below the boimdaiy symmetrical to 
point 3. It follows from the boundary conditions (26) and 
(38), both of which must be satisfied at the pomt 0, that the 
first derivatives in the coordinate directions are both zero 
at that point. On the basis of this fact, if the function 
yf''(v,0) is expanded in a two-dimensional Taylor’s series 
about point 0, the following finite-difference relations for the 
second derivatives are easily obtained: 

AV'„|o=4f'i-.| yl/z~l yl/o 
AV'm|o=2^'3-2v(''o-A"AV'„Io 

Here the symmetry property about the horizontal axis has 
been used to equate rj/'i to ^’3. Substitution of these relations 
into equation ( 20 ) for t 7=— 1 leads to the following finite- 
difference equation for the point 0: 

4(l-2F),A'i-| (1-2^] l^'o=0 (A3) 

FINITE-niFFERENCE EQUATIONS SPECIAL TO f'a 

The only finite-difference equation special to the problem 
for is the one used to terminate the field of computation 
at some vertical line on the left. As in the corresponding 
work for this equation is derived from an asymptotic 
solution of the boimdary-value problem valid for large 
negative values of 17. The derivation is parallel to that 
described in detail on page 16 of reference 2. 

The boimdary conditions which must be satisfied by yj/'s 
at large negative values of rj are shown in figm'e 16. The 
shaded section shows the anticipated variation of for 
constant tj. A solution of the differential equation which 
satisfies the given boundary conditions is 

a cos (g;)xV=^ (-217)”/*] 
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where K^z is the modified Bessel fvmction. of the second kind 
of order 1/3 and the On are appropriate constants. If only 
the leading term of this solution is used and the Bessel 
fimction is replaced by the first term of its asymptotic 
expansion, there results 


(A6) 

where F is the hypergeometric function and is 

defined by 

9 {B„-8y 

^ 8 rf 


Equation (A6) is suitable for use near the upper boundary 
(d^dayt^O). Near the sonic line (— ijeO, fei— =>) the 
following alternate form is available: 


22/3(l_f)l/6 




^ 3 3’ 1-f/^ 

(A7) 


If equation (A7) is evaluated on the sonic line, there results 

ql/3 

'FA.io, (as) 


rniv,e)=OcoB exp [-^ (-2,)^/^] 

As in the earher work, let A denote the lattice interval and 
/3 some lai’ge negative value of ij such that A//S<d. It then 
follows from the foregoing solution that, to a first order and 
for a given value of d, 

By substituting this relation into equation (Al), a finite- 
difference equation can be obtained which is valid for points 
on the line 17=— /3 and does not include any points to the 
left of this line (cf. eq. (39) of ref. 2). 

FINITE-DIFFEHENCB EQUATIONS SPECTAL TO 

The only equations special to the problem for 4'' a arise as 
a consequence of the condition along the upper boundary, 
where the values of 4 /j. are pjrescribed as a function of 17. 
Along most of the boundary, this condition can be met by 
substituting the prescribed values directly into finite-differ- 
ence equations of the type (Al) for points one interval below 
the boimdary. Because of the nature of the boundary values 
near 17=0, however, some change from previous procedures 
is necessary in the vicinity of the shoulder. Modification is 
also required in the equations used to terminate the field on 
the left. 

Points near shoulder of wedge. — ^From the known behavior 
of in the vicinity of the shoulder (see Guderley^s results, 
ref. 13, for the flow around a convex comer), it can be shown 
that the variation of along the upper boundary near 
17=0 must be of the form 

4''AiM=Mv,e^)=D{-vy/^ (A6) 

where Z? is a constant of proportionality. A singular solution 
of the differential equation (20) which is valid in the vicioity 
of the shoulder and which satisfies the boundaiy condition 
(A6) is also obtainable from Guderley’s results. This 
solution is, in the present notation. 


This result is in agreement with equation (32), which was 
developed from other considerations. It is apparent from 
equations (A5) and (A8) that a solution for will have a 
singularity in the first derivatives at the point 17=0, 6=6^,, 
Because of the foregoing singularity, a direct numerical 
calculation of i^'a might be expected to run into difficulty 
in the vicinity of the shoulder. Attempts along these lines 
led, in fact, to the unlikely result of negative lift over a 
small region of the profile just forward of the midchord. 
Reductions of the lattice interval to quite small values 
served merely to decrease the extent of this region. This is 
in contrast to the previous work for (and for as well), 
in which the singularity at the shoulder appears in the 
second derivatives. In that case, a sufficiently accurate 
solution for the unknown function could be obtained by 
direct calculation. In the present work, it was found 
necesdary to subtract out the singularity in the first deriva- 
tives according to the following procedure: , 

Let a function Sip ' a be defined such that 

. (A9) 

where a is a singular solution of the type given by equations 
(A6) and (A7). If the actual, numerically determined 
values of on the upper boundary are examined, it is 
found that for a small length of the boundary near the 
shoulder these values can be replaced to a good approximation 
by a )i-power variation of the form given by equation (A5). 
This is done, and the constant D is determined such that 
within this length of boundary \A'a,(i 7 , 5»)=^?''a(’ 7, O or 
sy_i=0. On this basis, a boundary-value problem for 
d\l/\ can be defined for a small region near the shoulder as 
shown in figure 17. The problem for within this region 
is solved jointly with the problem for in the remainder of 
the field. The two regions are fitted together by the use of 
overlapping lattices, much as is done in the case of a graded 
mesh (see ref. 15). The only difference is that equation (A9) 
must now be utilized to make the transition between the 
two lattices at all their common points. It is seen from 
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Figubb 17, — Boundary-value problem for 


figure 17 that conditions for 5^'^ on both the upper boundary 
and. sonic line are identical to the corresponding conditions 
for The finite-difference equations for the calculation of 
dtff'A can therefore be taken over directly from the previous 
■work. 

As nearly as one can judge -from experience •with various 
lattice spacings, results obtained by the foregoing process are 
quantitatively as well as quaUtatively rehahle. The primary 
source of error is in replacing the actual values of along 
the upper bmmdary by a }^-power variation. Since the re- 
gion over which this is done in the hodograph corresponds to 
a very small portion of the chord in the physical plane, 
errors from this source are probably small. 

Points far to the left. — The boundary conditions for 
at large negative values of 17 are shown in figure 18. BVom 

't'JivM V 



Figubb 18. — Boundary conditions on jP'a at large negative 


the asymptotic solution for the basic problem (eq. (38) of 
ref. 2 ), the expression for 4'' a along the upper boundary w 
foimd to be 

'yA(.v,0u)=l'eiv,0,r)=E{—ij)~^'* exp (— 

(AlO) 

where E]s & constant. 

Because of the nature of the boundary condition (AlO), it is 
not possible to write an asymptotic solution for for large 
negative 17 in a single term. For this reason, the procedure 
used previously to terminate the field of calculation at some 


location on the left cannot be apphed in the present case. 
An alternative procedure, more arbitrary in nature, can be 
de'yised by ‘writing in the form 

rAiv,e)=i''AivA)+srA{v,e) (All) 

where Stp'A is now defined by 5\j/' A{.‘nfi)='y A{i)fi)—4>' A{-nfiu) 
(see fig. 18). The attenuation of \{/'a in going from a point 
at 17 =— j3 to a point at 17 =— 73 — A is then foimd by treating 
each of the terms in equation (All) as an independent 
quantity. The attenuation of ^a(v,0u) is found from equa- 
tion (AlO) by a procedure similar to that used in obtaining 
equation (A4). The result is 

,A'^(-i3-A,d„)=[(l-A^ exp 

(A 12 ) 

To obtain a corresponding equation for 8^P'a, it is assumed 
that for a given value of 6 this quantity attenuates in the 
same manner as was pre'yiously found for One thus 
has from equation (A4) 

5-A'^(-^-A,0»)=[(l-^) exp 

(A13) 

Substitution of these expressions in equation (All) for 
17 =— 73 — A gives finally 

{[exp rA-w+rA(-M] (ai4) 

Since ^'^(— 73,^w) is a known quantity for any given value of 
73 , this equation can be used to terminate the field of calcula- 
tion in the same manner as was done 'with equation (A 4 ). 
The considerable element of arbitrariness in the derivation 
of equation (A14) can be tolerated since the over-aU solution 
is insensitive to changes in the left-hand portion of the field. 

SOLUTION OF FINITE-DIFFEKENCE EQUATIONS 

The techniques used to obtain a solution of the finite- 
difference equations for ^'a and yp'a were the same as those 
described in reference 2 for the basic solution In general, 
the graded lattice used for (see fig. 22 of ref. 2 ) was suit- 
able for the solution of For yj/'A, however, different 
gradations were necessary -with the smallest lattice spacing 
being used near the shoulder (point B of fig. 6 ). The value 
of i//'b at the intersection of the shock polar and the sonic 
line was chosen as 10,000 so that the pre'vdously obtained 
values of f could be used to proidde the initial guess for yl/g. 

In the course of the present work, a useful technique Avas 
found for locating regions of relatively large error in the 
numerical solution. By use of one form of Green’s theorem 
plus the differential equation ( 20 ), it can be shoAvn that 
around any contour enclosing a region in which equation ( 20 ) 
is satisfied the foUo-wing relation must hold: 

^ {2Tpl/'idr]-{-\l/'^)=0 (A16)' 
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In a numerical solution the line integral in equation (Al5) 
will not, except by rare coincidence, be precisely zero around 
any given contour. The amount by which it differs from 
zero may be taken as a rough measure of the adequacy of 
the numerical solution over the region within the contour. 
If the entire field of calculation is subdivided into a number 
of contiguous regions, it is thus possible, by evaluating the 


integral around each of the enclosing contours, to locate re- 
gions within which the error is relatively high. The solution 
in these regions can then be improved by advancing locally 
to a finer mesh. This technique was found to be of great 
help in the present work. It would probably be useful in 
other eUiptic boundary-value problems for wMch a relation 
analogous to equation (A15) can be obtained. 


APPENDIX B 


CALCULATION OF PLOW OVER REAR WEDGE IN PHYSICAL PLANE 


The procedure used to calculate the fiow over the rear 
wedge has been outlmed in the section METHOD OF 
SOLUTION. The fundamental operation is to determine, 
by stepwise methods, the initial rate of movement of the 
Icno-wn intersection points in the basic characteristics net. 
The methods which are used depend on the fact that these 
points are, by virtue of the basic characteristics construc- 
tion, points of fixed rj,6 (cf. eqs. (69) and (71) of ref. 2). 

The first step is to determine the initial rate of movement 
of those points at which the Mach lines of the basic charac- 
teristics net meet on the sonic line. For this purpose, con- 
sider equations (48), which give the initial rate of movement 
of a general point of fixed t],8. If these equations are spe- 
cialized to apply to points on the sonic line, the following 
relations are obtained: 

X'iO, 6)=-^ (Bla). 

F'(0,5)=^^^ _ (Bib) 

To wi’ite equation (Bla) the path of integration in equation 
(48a) is taken along the upper boundary from O to B (see 
fig. 5) and thence downward along the sonic line. The con- 
tribution of the portion from O to B is zero by virtue of con- 
dition (37). In applying these equations, the value of /„ is 
Imown from the basic solution. The integral in equation 
(Bla) is evaluated by mechanical integi'ation of a curve of 
munericaUy determined derivatives. Proper allowance is 
made for the singularity at the shoulder by integrating the 
singular solution analytically. The component rates of 
movement of the sonic point at the shoulder are both seen 
to be zero. 

The next step in the solution is to calculate the rate of 
movement of intersection points downstream of the sonic 
line. This is done by proceeding stepwise along consecu- 
tive downgoing characteristics. 

Consider three typical net points as shown m figure 19 (cf. 
also fig. 27 of ref. 2). The dashed lines represent the origi- 
nal position of the Mach lines through points a, b, and c, and 
the sohd lines represent their displaced positions correspond- 
ing to a small, finite angle of attack a. Since the intersection 
points in the Mach net are points of fixed r],d, the components 
of their displacement are given by aX' and ctF'. The slope 
of each segment of Mach line is taken, m accord with the 
procedures of reference 2, as the average of the slopes calcu- 
lated at the two end points. The slope calculated at each 
end point depends, in turn, only on the value of rj at that 
point (cf. eq. (68) of ref. 2). 



FiGtTEB 19. — Typical points in characteristics net. 


It is desired now to determine X' and Y' at point c in 
terms of X' and Y' at points a and b. Since the value of tj 
at a given net point is the same in the displaced and undis- 
placed positions, it follows from what has been said above 
that each segment of Mach line must retain its original slope 
after displacement. If this slope is denoted by m, the follow- 
ing relations are then readily obtained: 


_ Y'g—Y'b+ tutkX'i,— mgcX'a 




(B2a) 

(B2b) 


With these relations, it is a simple matter to calculate the 
initial rates of movement of successive net points on consecu- 
tive downgomg characteristics. For the first characteristic to 
be considered, point b is taken at the shoulder of the profile, 
where X' and Y' are both zero. Thus, X' c and Y' c for net 
points on this characteristic can be determined solely in 
terms of X'a and Y'a and the slopes riiac and rtitc- For the 
remainder of the downgoing characteristics, X\ and Y't are 
known from calculations along the characteristic immediately 
preceding. The actual calculations can be carried out in 
straightforward tabular form. 

The foregoing procedure enables the calculation of X' and 
Y' for all net points except the ones originally at the surface 



566 


REPOKT 1180 — ^NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


of the rear -wedge. For these points, consideration must be 
given to the required boimdary condition at the surface. 
This boimdary condition is 

0(X,+O;a)= — (Offi+a) (B3) 

from which it follows that 

0'(Z,+O)=-1 (B4) 


The problem now is to determine X' and T' at the surface 
of the wedge in such a way that equation (B4) is satisfied. 
To do this equation (53b) is first specialized to the surface 
of the wedge, where it is readily shown that Xe = Y^=Q. In 
view of condition (B4), there results 

X'(.v,—&w)=Y0(.‘n,—6m) (B5) 

The value of Y' at points originally on the surface of the 
wedge is thus fixed directly by the basic solution. The cor- 
responding value of X' can be found from a construction 
analogous to that of figure 19 and is 


X' = 


Y\—Y'c—m.aoX'a 


—m.. 


(B6) 


The point c is now the point originally on the surface of the 
wedge 6- ®-j is as given by equation (B5)), and the re- 
maining notation is the same as in figure 19. 


Application of equation (B5) requires the knowledge of 
Yii^, — O, which in the case of the wedge profile is equal 
to 1/tfy. Evaluation of the latter derivative can be carried 
out directly from the basic Mach net, but the procedures are 
cumbersome and inaccurate. A better method is to use the 
equations of motion (cf. eq. (6) of ref. 3) to express 6y in 
terms of ijx- Following this procedure, one obtains finally 


Y$(,r],—e„)= 


{2d^yi^ 

2vi^,0) vx(^,+0) 


(B7) 


The quantities Ij and Ijx which appear here are easily eval- 
uated from the basic solution for the chordwise distribution 
of^. 

The preceding equations enable the calculation of the initial 
rate of movement X' for points originally on the surface of 
the rear wedge. The final step is to determine the cor- 
responding distribution of lift. For this purpose, equation 
(53a) is specialized to poiuts on the rear wedge to obtain 

(A,_+0) = -X' (^,- O 

which, in view of the Jboundary conditions, can be sho^vn to 
be equivalent to 

ij'(X,-t-0)= — X'{r), — 0») i?x(X,-l-0) (B8) 


The distribution of hft is then obtaiued from equation (57). 


APPENDIX C 


SOLUTION OF PROBLEM FOR COMPLETELY SUPERSONIC FLOW 


CALCULATION OF UFT-CURVB SLOPE AND 
CENTER OF UFT 

If conditions are such that 2'/®= 1.260 (corresponding 
to 1; cf. eq. (13)), then the basic flow over the profile at 
zero angle of attack is completely supersonic. The solution 
for the lift-curve slope and center of lift at a vanishingly small 
angle of attack can then be carried out analytically as follows; 



Consider a completely supersonic flow about the double- 
wedge profile at a small angle of attack. In the physical 
plane the flow field has the weU-kno-wn appearance sho-\vn on 
the left in figure 20. The corresponding hodograph of the 
flow along the upper surface, in terms of the normalized 
small-disturbance variables tj and 6, is sho-wn on the right. 
The quantities and a are, as before, the half-angle of the 


e 



\ 


Figxteb 20.— Conditions on airfoil in completely supersonic flow. 
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wedge and the angle of attack (also normalized). Except 
for a small range of just above 1.260 (see below), flow 
conditions must be constant along each of the segments 1 
and 2 of the upper surface. In the hodograph each of these 
segments is thus represented by a single point located as 
shown. It is apparent that for a given value of d„, the 
speeds and which are the primary unknowns in the 
problem, are fimctions solely of the angle of attack a. 

To find the lift-curve slope and center of lift it is necessary 
first to find the derivatives ri'i=(,dtiiJdat)a-o and ij't= 
idi) 2 lda)a^{j. This can be done with the aid of the equations 

for the shock polar 

0=(1— )7)Vi+^ (Cl) 

and for the downgoing characteristic 


23/3 

0=constant — ^ 
o 


( 02 ) 


To find rj'i, one must utilize the boundary condition 
6i=0a—a. Substitution of this condition into equation 
(Cl) provides the following implicit equation for iji: 


dy,— a= (1 — iiO 

Differentiation of this equation gives 


(C3) 


da l~4~3i;i 
dvi 2-\Zl-f-)?i 


From this it follows that 


V 1 = 


2Vl4~?i 

1 + 3 i 7 i 


(C4) 


where, as in the main text, the bars denote the value of 
7 ]i at a=0. The value of iji can be foxmd in terms of the 
parameter da by solving equation (03) for i/i with a set 
equal to zero. The result, obtained through standard 
methods for the solution of cubic equations, is 


Vi=- 




2 cos 


■7T — ip 


■1 


(C5) 


whore 


v?=axc cos 


(I 

\4 V 2 


To find rj'i, equation (C2) for the downgoing character- 
istic is first specialized so as to pass through the point 1. 
This gives 

03/2 

e=(eg-a)+^ 

Substitution of the boundary condition 63 —— 6 a—a then 
provides the result that 




II Compare eqoatlon (67) of relorenco 2. 
308CC5— 00 3T 


Takihg the derivative with respect to a, one obtains finally 
at Q==0 

K - v/3 (C6) 

■ 




where i]'i is given by equation (C4).and by equation (C5). 

Since the value of 77' is constant on each segment of the 
profile, the lift-curve slope is easily found from equation 
(57) and is 

[(T-fl)MJ(i/c)]^^(^)^_^=2(20„)>«(/i+7,',) 
Substitution from equation (C6) gives 
[(7-K)MJ(i/c)]>/®(^)^_^=2(2fl„)>«77' 


r X 1/2 “I 

11 

1 / q \i/3 


(C7) 

The moment-curve slope, for moments taken about the 
leading edge, is found to be 

(2«»)‘'^(’7'f+37;'2) 


or 






1 + 


3? 


«.i/s 


( 


» \l/3 


3 

V2' 


(C8) 


The position of the center of lift is given accordingly by 

1 I 3^1^/'“ 


0 ), 


{dci/do^a~~o 4 ^ 






m 


[^i^'*+(3/V2)0„]*« 


In equations (07) and (08), the first term inside the brackets 
represents the contribution of the front wedge, the second 
term that of the rear. 

Equations (07) and (09) are the basis for the curves 
shown in figures 11 and 13 for values of 1.260. The 
results show certain curious features when the flow over the 
front wedge is just sonic (i. e., ^i=0, 0 h,= 1, f„=1.260). 
These are as follows: 

(a) The lift contributed by the rear wedge is zero (see 
eq. (07)). 

(b) The center of lift is at the quarter-chord point (follows 
from statement (a) plus the condition of uniform lift on the 
front wedge; see also eq. (09)). 

(c) The rate of change with respect to is infinite both 
for the lift-curve slope of the complete profile and for the 
position of the center of lift (follows from differentiation of 

'll ni .1 

(SO) m-.: . 

d ihm'i 
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eqs. (07) and (09)). These results ai-e associated in every 
case with the behavior of the lift calculated for the rear 
wedge. 

ESTIMATION OF LOWER LIMIT FOR CONSTANT SPEED 
ALONG REAR WEDGE 


tan X 


(2d^y'^ 
Vl + ’7l 


Au analogous expression for the Mach angle /x is given by 
equation (68) of reference 2 and is 


The features just enumerated, though having a certain 
curiosity in themselves, cannot be accepted as completely 
correct. Because of interaction effects between the shock 
wave from the bow and the expansion fan from the shoulder, 
the fundamental condition of constant- speed at the surface 
of the profile wiU not be satisfied alqng the rear wedge until 
the value of is somewhat greater than 1.260. Until then, 
disturbances reflected from the shock wave wfll reach the 
rear wedge and cause a slight decrease in speed toward the 
trailing edge. This effect will cease when the forward most 
reflected Mach wave just touches the trailing edge. The 
exact value of at which this condition wiU be met is 
difficult to determine. An upper bound can, however, be 
estimated as follows: 

Consider the basic flow field (a=0) over the upper half of 
the profile when the first reflected Mach wave just strikes 
the trailing edge. Figure 21 shows such a flow field as it 
would appear in transonic similarity form (cf. pp. 5-6 of ref. 
2). In drawing the figure a special assumption has been 



Figure 21. — Macli-liiie pattern assumed to esist when first reflected 
Mach wave strikes trailing edge. 

introduced beyond those implicit in the small-disturbance 
theory; namely, that the first reflected Mach wave is straight 
and has an angle of inclination /t equal to that of thq first 
wave in the expansion fan. With this assumption, the cor- 
responding value of is easily determined. Since the 
reflected wave must actually be curved downstream, the 
value so determined will be greater than the correct value 
for the required condition. 

On the basis of figure 21, the following equation can be 
written between the shock angle X and the Mach angle jx: 

tan X=itan ft (ClO) 

A relation between the shock angle X and the speed 
region behind the shock can be obtained from equatioj(i^’05&) 
and the known properties of the shock polar. The 


tan 


(20„)i« 


Substitution of these relations into equation (OlO) and 
solution for gives 



The accompanying value of foimd from equation (03) 
with o;=0, is 

0„=O.9685 

This corresponds, according to equation (13), to 

f„=1.287 (Oil) 

Thus, for values of between 1.260 and some limit less 
than 1.287, the results of equations (07), (08), and (09) are 
not exact insofar as the contribution of the rear wedge is 
concerned. It can be reasoned that in tins range an exact 
solution would indicate more lift for the rear wedge than 
does the present analysis. 
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